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Abstract 

We employ effective interaction potentials to study the equilibrium structure and phase behav- 
ior of highly asymmetric mixtures of star polymers. We consider in particular the influence of the 
addition of a component with a small number of arms and a small size on a concentrated solution 
of large stars with a high functionality. By employing liquid integral equation theories we examine 
the evolution of the correlation functions of the big stars upon addition of the small ones, finding 
a loss of structure that can be attributed to a weakening of the repulsions between the large stars 
due to the presence of the small ones. We analyze this phenomenon be means of a generalized 
depletion mechanism which is supported by computer simulations. By applying thermodynamic 
perturbation theory we draw the phase diagram of the asymmetric mixture, finding that the ad- 
dition of small stars melts the crystal formed by the big ones. A systematic comparison between 
the two- and effective one-component descriptions of the mixture that corroborates the reliability 
of the generalized depletion picture is also carried out. 

PACS numbers: 61.20.-p, 61.20.Gy, 64.70.-p 
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I. INTRODUCTION 



Mixtures whose constituent particles show a high asymmetry in sizes are quite common in 
soft matter physics. As a matter of fact, all soft-matter systems are at least two-component 
mixtures, as they are typically suspensions or dispersions of mesoscopically-sized colloidal 
particles in a microscopic solvent. For many practical purposes, though, it suffices to model 
the solvent as a continuous medium and then an effective, one-component description of the 
suspended colloidal particles is sufficient. The phenomenology is much richer when more 
than one species of colloids is dispersed in the solvent and also there the asymmetry in the 
sizes of the two kinds of colloidal particles can be much higher than the one encountered in 
atomic or molecular fluids. In the recent past, a great deal of attention has been paid to the 
investigation of model colloid-polymer mixtures, in which the two species are hard colloidal 
spheres and soft, flexible polymer chains The bulk of the theoretical analysis of such 
systems is carried out within the framework of the Asakura-Oosawa (AO) model 0,0,0], in 
which the polymers are modeled as ideal, interpenetrating spheres that experience a hard 
repulsion towards the colloids. Another popular system that has attracted a lot of attention 
recently are binary hard-sphere mixtures of various size ratios 00, zl si in both of 
those cases, attention is usually focussed on the influence of the smaller component on the 
structural and phase behavior of the larger one Q]. Demixing phase transitions and their 
competition to the crystallization transition of the large hard spheres have been an issue 
of intensive investigations in the past QHHflfiyHQ with cu rr e„t research 
stee ring in the direction of the study of interfacial and wetting properties of such mixtures 
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A convenient concept that has helped shed light into the phenomenology of such asym- 
metric mixtures is that of the effective, depletion interaction between the hard spheres, which 



Hletic 
Q. 



is mediated by the smaller component (2J, 35]. In the case of the AO-mixture, the depletion 
interaction is purely attractive and has the range of the size of the added polymer. For 
binary hard-sphere mixtures, the effective depletion potential displays oscillatory behavior 



due to correlation effects 



36 



371 ]. Interpolating between the two extremes of ideal and hard 



additives are star polymers of varying functionality, whose depleting effects on hard spheres 
have been investigated both theoretically QQ and experimentally Q. 
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The notion of depletion is almost exclusively invoked whenever the large particles are 
hard colloids. Nevertheless, it can be expanded in its interpretation to account for the 
modification of the properties of the large particles in the presence of smaller ones also for 
arbitrary kinds of interactions between the constituent particles. There is relatively little 
done in this direction, however, with the exception of the derivation of effective potentials 
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in Yukawa mixtures [4l[ and in mixtures of star polymers and linear chains 
last case, it has been shown that the depletion mechanism of the chains on the stars can 
account for the experimentally observed melting of the star-polymer gel upon addition of 
linear polymer. In this paper, we turn our attention to two-component mixtures in which 
all particle species interact by means of soft potentials and, in particular, to mixtures of 
two kinds of star polymers: large ones with a high number of arms and small ones will a 
low arm number. All species interact via logarithmic- Yukawa pair potentials. We find that 
in this case the depletion mechanism of the small stars on the big ones has the effect of 
reducing the repulsive potential between the latter and thus it brings about a melting of 
the colloidal crystal formed by the large stars. Concomitant to this effect is a partial loss 
of correlations between the centers of the big stars, manifested in a drastic lowering of the 
peak height of their partial structure factor. Upon addition of a sufficiently large quantity 
of depleting agents, even an effective attraction between the large stars shows up, resulting 
in a demixing spinodal between the two species. 

The rest of the paper is organized as follows: In Sec. |H]we present the pair potentials and 
the full, two-component description of the mixture, examining the effects of the depletants 
on the structural correlations of the big stars. In Sec. we formally trace out the small 
stars and examine the resulting effective, one-component interactions between the big ones. 
This effective potential is employed, in turn, in order to draw the phase diagram of the 
system in Sec. IIV[ where thermodynamic perturbation theory is used for the calculation of 
the Helmholtz free energies of the fluid and solid phases. In Sec. E]we carry out a comparison 
between the one- and full two-component descriptions of the mixture and demonstrate the 
validity of the former, whereas in Sec. IVII we summarize and draw our conclusions. 
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II. TWO-COMPONENT DESCRIPTION 



We consider binary mixtures of star-polymers which differ in terms of their sizes and 
arm numbers (functionalities). The system consists of N\ stars of corona diameter a\ and 
functionality j\ and N 2 stars, characterized by o 2 and f 2 , in a volume V. We first calculate 
the properties of the binary fluid. To obtain information about the pair correlations between 
the constituent particles, we describe the system using the full two-component picture for 
the mixture of the two different star-polymer species. The structural quantities we calculate 
are used as input for the mapping onto an effective one-component system in Sec. IIHI We 
define the size ratio of the different species as q = cf 2 /(J\ < 1. Let pi = N^/V {i = 1, 2) be 
the partial number densities of the two species. 

We start from the effective pair potentials between the mesoscopic particles, having traced 
out the monomer and solvent degrees of freedom. The effective interaction between the star- 
polymers diverges logarithmically with the center-to-center distance r as r — > 0, as derived 
by Witten and Pincus A full expression for identical star-polymers, which is valid for all 
star separations, has been derived and verified by neutron scattering and monomer resolved 
molecular simulation Q| • The pair potential is given by an ultrasoft interaction which 
shows logarithmic behavior for small distances and an exponential Yukawa-type decay at 



large star-star separation 
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In the case of mixtures we need an expression for the effective interaction between star 

n 

polymers in an athermal solvent that differ in their sizes a 2 [70] and functionalities f±, 
/ 2 , as a function of their center-to-center separation r. In this work we use the effective pair 
potential which was put forward by means of field-theoretical arguments and confirmed by 
molecular dynamics computer simulations in Ref. [47], namely: 
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Moreover, f3 = (ksT) 1 is the inverse temperature, with &b being Boltzmann's constant. 
Since all three interactions are purely entropic, the j3Vij(r) are independent of the tempera- 



ture. For i = j the potential reduces to the interaction of identical star-polymers which was 
introduced in j^J . In what follows, we fix the functionality of the large stars to f\ = 263 in 
order to make contact with recently performed experiments [42] in which smaller polymeric 
entities were used as additives in gelated solutions of the large stars in order to examine 
their overall influence on the rheology of the mixture. This functionality is large enough for 
ihe star polymers to crystallize into a fcc-structure roughly at their overlap concentration 
4||. For the small stars, we considered functionalities f2 = 16 and 32 and size ratios q in 
the range between 0.1 and 0.3. 

The pair structure of the mixture can now be calculated using the Ornstein-Zernike (OZ) 
equations for binary mixtures together with the two-component Rogers- Young (RY) closure. 
The pair correlations of the system are described by three independent total correlation 
functions hij(r), i < j = 1,2, since the symmetry with respect to the exchange of indices 
dictates hij(r) = hji(r) for % ^ j. In addition, we have the same number of direct correlation 
functions Qj(r). The Fourier transforms of these quantities are denoted by hij(k) and (%(&), 
respectively. 

For multicomponent mixtures, the OZ relation takes the form 0, |3] 

H(k) = C{k) + C{k)-B-H{k), (3) 
where H(&) and C(k) are symmetric v x v matrices with 

[H.(k)}ij = hij(k) and [C(k)]ij = (%(&). (4) 
D is a diagonal v x v matrix with 

\D]ij = PiSij. (5) 

From Eq. (J3J) we obtain three independent equations for the six unknown functions hij(k) 
and Cij(k), i,j = 1, 2. In order to obtain a solvable system, we need three additional closure 
equations between these functions. The Rogers- Young closure for multicomponent mixtures 
reads as [34] 



gij(r) = exp[-/3V^(r)] 
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where gij(r) = hij(r) + 1, Jij{r) = hij(r) — Cfj(r) and V^-(r) is the pair interaction between 
particles of species i and j. The "mixing function" fij(r) is defined as 

fij( r ) = 1 -exp(-ayr). (7) 
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FIG. 1: The partial (a) radial distribution functions and (b) static structure factors for species 1 
(big stars) in a mixture with small ones. The density of species 1 is p\o\ = 0.05. The structure 
has been calculated using the two component OZ equations with the RY closure. The plotted lines 
are for different densities p2, as indicated in the legend. In this case, q = 0.1 and fi = 32. The 
partial structure factor grows for k — > with increasing p2, as we approach the demixing spinodal 
line of the system. 

Usually, the same self-consistency parameter a = ctij is used for all components of the 
mixture. This allows to fulfill one thermodynamic consistency requirement, namely the 
equality between the "virial" and "fluctuation" total compressibilities of the mixture. Multi- 
parameter versions have also been proposed |6j , invoking thermodynamic consistency for the 
partial compressibilities of each species. For a — > Eq. (jUJ) reduces to the Percus-Yevick 
(PY) and for a — ► oo to the hypernetted chain (HNC) multicomponent closures. When 
dealing with star-polymers, which feature a soft repulsion of relatively short range, neither 
the PY nor the HNC closure are adequate to capture the details of the correlation functions 
with high accuracy, therefore employing the full RY closure is essential 49]. 

In our work we solve the OZ equation with the RY-closure [Eqs. ©-(0)] for the two- 
component mixture. The effective interactions between the star polymers are given by Eq. 

The thermodynamic consistency of the RY closure was obtained by using a single 
parameter a. The structure of the binary mixture can be described either by the partial 
radial distribution functions gij(r) = hij(r) + 1 in real space or by the three partial static 
structure factors Sij(k) = 5ij + y/PiPjhij(k) in wavenumber space. The structure factors 
are relevant in comparing with experiments, because they can be measured via scattering 
techniques. 

In Fig. ^a) we show results for the radial distribution function gu(r) between the large 
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FIG. 2: Same as Fig. ^ but now for big star density pi<jf = 0.3. By increasing the density of the 
smaller component, the structure of the fluid diminishes. 



stars in a dilute solution, and its evolution upon increasing the concentration of small addi- 
tives with / 2 = 32 and q = 0.1. Although for very small concentration of smaller stars the 
function gu(r) has a relatively structureless shape, it rapidly develops a pronounced peak 
when p 2 is further increased. This is a first indication of clustering of big stars, which has 
its physical origin in some effective attraction induced by the small component. One phys- 
ically expects that when this attraction becomes sufficiently strong, a demixing transition 
between the two species will take place. This hypothesis is corroborated by the evolution 
of the corresponding structure factor Sn(k), shown in Fig. IHb). A fluid-fluid demixing 
binodal is indicated by the divergence of all partial structure factors in the long-wavelength 
limit k — > 0. As can be seen in Fig.^b), a growth of the k — > 0-limit occurs upon increasing 
p 2 . The existence of a demixing spinodal will be confirmed in Sec. IIVI where we draw the 
phase diagram of the mixture. 

We now examine the effect of the additives at the complementary regime of high concen- 
tration of large stars, and in particular slightly above their overlap concentration, in which 
the latter are in a thermodynamically stable crystalline state j^l or in a dynamically ar- 
rested gel state [50]. We derive the partial structure factor Sn(k) of the (metastable) fluid 
in the absence of small stars and monitor its evolution as p 2 is increased. Representative 
results are shown in Fig. El In Fig. Efa), already a loss of correlations in real space can be 
discerned, as witnessed by the broadening and lowering of the coordination peaks in gn(r). 
Moreover, the large stars approach closer to each other upon an increase of p 2 , an effect 
that can be interpreted as a weakening of the strength of their mutual repulsion. As can be 
seen in Fig. |2fb), the principle peak height of the structure factor of species 1 diminishes as 
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the density of the smaller component is increased. The Hansen- Verlet criterion |5l( states 
that a fluid solidifies when the maximum of the structure factor exceeds the threshold value 
Sth(^max) = 2.85. Therefore, the diminishing of structure in the system is a first indication 
for the melting of the crystal of big star-polymers by addition of the smaller species. In Fig. 
I2fb) it can be seen that the first peak of the structure factor is bigger than 3 for p2®\ = 0. 
Already small densities of the smaller component lead to a drastic decline of the peak height, 
a finding that is in line with recent experimental and theoretical results on mixtures of star 



polymers with linear chains 
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III. EFFECTIVE ONE-COMPONENT DESCRIPTION 

We now wish to put the assumptions regarding the influence of the additives on the 
effective interaction of the big stars into a concrete test, by calculating an effective potential 
V e g (r) between the latter in the presence of the former. To this end, we carry out the mapping 
of the two-component system onto an effective one-component description, in which the 
degrees of freedom of the smaller star-polymers have been traced out. The interactions cause 
spatial correlations of the density of small stars in the vicinity of two big ones, influencing 
thereby the shape of the resulting generalized depletion interaction. There are different 
methods to obtain these effective interactions. All of them omit many-body forces and 
reduce the interaction to an additive pair potential, we will confirm however that many- 
body effects only play a minor role. 

Instead of the so-called 'system representation', in which the two densities p\ and p 2 in 
the mixture are given, we now switch into the more convenient 'reservoir representation' 
Pi and p\. Since the effective interaction between the large stars depends rather on the 
chemical potential Uo of the small ones rather than on their density P2, this description 
is more convenient The reservoir is a system consisting of pure small stars and their 
density there, p r 2 , is determined by the requirement that the partial chemical potentials p 2 
in the real system and p 2 in the reservoir are equal. Clearly, due to the finite value of the 
density p\ in the system, it must hold p 2 ^ p 2 - The mapping between the two densities, 
depending parametrically on the big star density pi, will be carried out in Sec. |V| 
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FIG. 3: Comparison of Monte Carlo simulation (MC), inversion of Ornstein-Zernike equation (OZ), 
and superposition approximation (SA) results for the depletion forces between the big stars, (a) 
Results for q = 0.3, = 32 and p\a\ = 0.027; (b) q and f 2 same as in (a) but p\a\ = 0.081; (c) 
for q = 0.2, f 2 = 32 and p\o\ = 0.08. 

A. Simulations 



The most accurate way to obtain the effective ^interactions between the big star-polymers 
is a computer simulation of the mixture 
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We place two big stars of species 1 at 



fixed positions Ri and R2 along the diagonal of the simulation cube, so that their common 
center coincides with the center of the cube. They are surrounded by the smaller species 
that move according to the forces dictated by the effective interactions of Eq. Since we 
have only two big star-polymers in our simulation box, the density is p\ —>■ 0. Therefore 
the simulation provides directly the sought-for effective force as a function of the reservoir 
density p\. 

We use standard NVT Monte Carlo simulation with periodic boundary conditions and 
minimum image convention. The length of the cubic simulation box is L = 5<Ti, so that the 
number of small stars in the simulation results which are shown in Fig. El is between 125 
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and 1250. For each particle up to 5 million Monte Carlo steps are calculated, where the 
maximum displacement of the particles is chosen in such a way that half the steps will be 
accept. The force is then measured after every 1000 simulation steps. Due to the presence of 
the second big star, the density distribution around each star is not spherically symmetric. 
This leads to an average nonvanishing force between them, which is mediated by the small 
star-polymers and parallel to the vector Ri 2 = R2 — Ri- Due to the symmetry of the system, 
the components perpendicular to R12 have to vanish. The force Fi acting on the particle at 
Rx can then be calculated by averaging over the simulation results, namely: 



where rj are the star-polymer positions of species 2 and (• • • ) R denotes the statistical 
average, taken under the constraint of constant R12. Clearly, the effective force satisfies 
the relation F 2 (i?i 2 ) = — Fi(i?i 2 ). We further define the depletion force Fd ep (-Ri2) as the 
difference between F 2 (i?i 2 ) and the direct force F dir (i?i 2 ) between the two stars due to their 
direct interaction potential V n (R 12 ). The magnitude of the depletion force F dep (i2 12 ) is then 
given by 



Accordingly, the total effective interaction between the big star-polymers in a sea of the 
smaller species is the sum of their interaction potential Vu(r) and the depletion potential 



A large number of long simulation runs is required to to obtain accurate depletion forces 
with good statistics, which renders this approach inefficient if ones needs to calculate V e s(r) 
for arbitrary values of q, f 2 and p\. Thus, we resorted to approximative theoretical methods 
to calculate the effective interactions and used the simulation results at selected parameter 
combinations in order to put the theoretical approximations into test. The two theoretical 
approaches invoked in this work are the inversion of the Ornstein-Zernike equation and the 
superposition approximation, which are presented below. 






(9) 




V eS (r) = V 11 (r) + V dcp (r). 



(10) 
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B. Inversion of the Ornstein-Zernike equation 



The effective potential can be obtained by inversion of the two-component OZ equation 



results in the limit of low density of big stars 
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531 ] . It can be shown from diagrammatic 



expansions in the theory of liquids [57[ that the radial distribution function g(r) of any fluid 
whose constituent particles interact via the pair potential V(r), reduces to the Boltzmann 
factor g(r) = exp[—/3V(r)] in the low-density limit. The effective interaction between the 
big stars depends on the reservoir density p\ of the smaller component. The interaction 
can be obtained by solving the full two-component OZ equations with the RY closure for 
different small-component densities p 2 in the limit p\ — > 0; due to the latter limit, it then 
also holds p 2 = p\. The radial distribution function gu(r) can then be inverted to yield the 
effective potential as 

PV eS (r) = -\im\n[ gil (r; Pl ,pl)]. (11) 

Thereafter, the depletion force Fd ep (r) can be calculated as Fd ep (r) = — d\y e g(r) — Vn(r)]/dr 
and compared them to the simulation results of the preceding subsection. Selected compar- 
isons are shown in Fig. |21 where it can be seen that the inversion of the OZ relation yields 
very reliable results. We emphasize here that the approximate character of the OZ-inversion 
technique lies exclusively in the approximations involved in solving the two-component in- 
tegral equation theories, i.e., in the Rogers- Young (or any other chosen) closure relation. 
Otherwise, the method is based on the exact statement that the radial distribution function 
of a one-component system at low densities is equal to the Boltzmann factor of the associated 
pair potential and hence the agreement of the inversion method with the simulation results 
comes as no surprise. It rather corroborates the fact that the two-component RY-closure is 
very accurate whenever one deals with soft, repulsive interactions, a result already seen in 
the case of mixtures between hard spheres and star polymers 



C. Superposition approximation 

Another possibility to derive the effective interaction is the superposition approximation 
(SA) 5^|. If the density distribution p 2 (r;R!,R 2 ) of the small stars around two big stars 
held fixed at positions Ri and R 2 is known, then the depletion force in the low-density 
limit can be calculated by a simple integration. The density p 2 (r;Ri,R 2 ) is proportional 
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to three-body distribution function gn 2 (Ri, R2, r), which expresses the probability density 
of finding a particle of species 2 at position r, given that two particles of species one are 
fixed at positions Ri and R 2 . This function is in general unknown; a usual procedure is to 
approximate it by the product of pair distribution functions (5^. 

We consider two big stars at the positions Ri and R 2 and choose, without loss of gen- 
erality, Ri = 0. Let the distance between the particles be Ru- The surrounding smaller 
star-polymers have the density p 2 (r; Ri, R 2 ). By taking the average for fixed Ru we obtain 
the depletion force as 



Fdep(Rl2j 



2n I r 



1 



,dV l2 (r 



-dr 



dr 

p 2 (r; Ri, R 2 )a>da;, 



(12) 



where u = cos 9. 



P 2 (r;R,Jtj) 




FIG. 4: A sketch of two big stars at a distance R\2 with p 2 (r; Ri, R 2 ) denoting the density of the 
smaller stars at r. The density distribution depends on the positions of the two big star-polymers. 

Since p 2 (r;R 1; R 2 ) is in general not known, at this point the exact density distribution 
has to be approximated. The density distribution around two big stars is replaced by the 
product of the density distributions around two isolated star-polymers at the positions R x 
and R 2 , respectively. The SA then reads as 

p 2 (r; R 1; R 2 ) « p 2 g 12 (\r - R 1 |)^ 12 (|R 2 - r|), (13) 

where p 2 is the reservoir density, again identical to the system density for the situation at 
hand, since only two big stars are considered in the thermodynamic limit and thus pi = 0. 
The functions gi 2 (|Ri — r|) are the radial distribution functions of small star polymers 
surrounding a single large one. Therefore, they can be obtained in the pi — > 0-limit of the 
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two- component OZ equations. Using simple geometrical considerations, we obtain |r — Ri 



The results we obtain by this method are compared to Monte Carlo and to those derived 
from inversion of the Ornstein-Zernike equation in Fig. EJ The results are very similar to 
the ones we obtain by inverting the OZ equation and both approximations yields reasonable 
agreement with the simulation data. Therefore, it is possible to choose the results of either 
approximation for calculating the phase diagrams and we expect that only minor quantitative 
differences will be seen by employing the one or the other theoretical approach. We have 
chosen to work with the effective interactions that result from the inversion of the OZ- 
equation since the latter is based on an exact statement, whereas the SA has an approximate 
nature. 

D. Effective interactions 

We have chosen to employ the effective interactions that result from the inversion of 
the OZ-equation since the latter is based on an exact statement, whereas the SA has an 
approximate nature. Representative results are shown in Fig. For the lowest size ratio, 
q = 0.1, we see that irrespective of the functionality of the depletants (f 2 = 16 or 32), the 
following scenario materializes: as p\ increases, first a weakening of the repulsions takes place, 
followed by the development of an attraction between the stars at sufficiently high reservoir 
densities, see Figs. E^a) andEfb). These findings provide a possible physical realization 
of the recently- prop osed model ultrasoft repulsion potentials that are accompanied by an 
attractive part 59]. This attraction is more pronounced for f 2 = 32 than for / 2 = 16, if 
one compares two systems with equal density p\a\. This result is not surprising, since the 
f 2 = 32-stars exert a higher osmotic pressure on the large ones than the f 2 = 16-stars and 
can therefore reduce the direct repulsions and induce attractions more efficiently. 

Novel features in the effective potential appear for higher size ratios, q = 0.2 and q = 0.3. 
As can be seen in Fig. Efc)-(f), an oscillatory structure appears in the effective potential, 



V -^12 + r2 — R\ 2 ruj. Finally we obtain for the depletion force in the SA the expression: 




(14) 



13 




FIG. 5: The effective potential V e s(r) between two big star polymers in the presence of a sea of 
small ones. The various combinations of parameters regarding the density, functionality and size 
of the additives are shown in the legends. Notice the development of a strong attractive part in 
the interaction for the case q = 0.1. 

which is akin to that seen for hard-sphere mixtures of two different sizes. Contrary to this 
case, however, a deep attraction between the big stars does not develop and, therefore, it 
seems that a demixing transition between the two species does not exist when the sizes of the 
two stars become more and more similar. In all cases, however, the range of the repulsion 
decreases due to the depletion effect, i.e., the big star polymers appear, in the presence of 
the small ones, to be softer than they are in a pure solvent. Another possible interpretation, 
to be elaborated on in what follows, is that the big stars appear to be 'smaller', i.e., they 
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FIG. 6: The two possible effective hard sphere diameters d and do pertaining to the effective 
interaction V e &(r) between big stars, against the reservoir density p\- The parameter combination 
here is fi = 32 and q = 0.1. As explained in the text, d is calculated using the full effective 
interaction and do only for the reference part. 

acquire a reduced effective hard sphere packing fraction as a result of the depletants. Since 
the star-polymers then need less space, they become more mobile so the solid can melt. This 
property will be discussed in more detail in Sec. IIVB1 

IV. PHASE DIAGRAMS 
A. Hard sphere mapping 

In order to trace out the phase diagram of the mixture in the (pi, /^-representation, we 
first perform a mapping of the effective one-component interaction V e g(r) between the big 
stars onto an effective hard-sphere system of diameter a. Clearly, the latter depends on 
the reservoir density p\ as well as on the system parameters q and For the purposes of 
performing the mapping in a physically meaningful way, we distinguish between two cases. 

First, we consider the case in which Kff( r ) is either free of attractive parts or positive 
definite or, at most, it contains negative parts not exceeding a small fraction of k^T in 
magnitude. In this case, it is physically meaningful to identify a with the Barker-Henderson 
hard sphere diameter of the full effective interaction V e ff(r), d, denned as \&\ : 



Most of the curves shown in Fig. fall into this category. An important exception are the 




(15) 
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curves pertaining to p\cr\ = 0.1 in Fig. Efa) and to P2 a 2 = 0-05 in Fig. 0(b) . For these com- 
binations, and also for all others at even higher reservoir densities, a deep negative minimum 
appears in V e s(r) and application of Eq. (JT5J) to such cases would lead to unphysically small 
and even negative effective hard sphere diameters. For such combinations, it is physically 
appealing to separate the effective potential V e s(r) into a purely repulsive part V (r) and 
a perturbation part V peTt (r), by truncating and shifting upwards the full interaction at the 
deepest minimum 6(|. In this second case, it is pertinent to define another effective hard 
sphere diameter, d , that is associated with Vq(t) only and is calculated again from the 
Barker- Henderson recipe, namely 

POO 

d Q = / {l-exp[-/5K,(r)]}dr. (16) 
Jo 

In attempting to choose and match between the two possible hard sphere diameters, d 
and do, we are confronted with a technical difficulty. The evolution of the potential V e g(r) 
with p\ is continuous and the appearance of negative minima is in general accompanied by 
a soft repulsive barrier after the minimum. The effective hard sphere diameter, on the other 
hand, has to be a continuous function of p 2 , so as to avoid unphysical jumps of the phase 
boundaries in the phase diagram. In Fig. |H]we show a typical result for the dependence of d 
and do on p\. For low values of p T 2 , where V e s (r) is purely repulsive, d is a meaningful measure 
of the effective hard-sphere diameter. On the other hand, at high values of p r 2 , where a deep 
attraction between the big stars effectively sets in, it is do that most realistically captures 
the physics of the repulsions. The two curves cross at some point and, in order to guarantee 
both the continuity of a as a function of p\ and its correct asymptotic behavior for small 
and large values of p 2 , we choose 

a = max{d, do}. (17) 

It is then clear from the discussion above that for the perturbation part, V^ ert (r), of the 
interaction, it holds 

V pevt (r) = iid>d . (18) 

The phase diagrams can now be calculated using standard first-order perturbation theory 
(^1 and taking Eqs. (fT7|) and ()18j) into account. We do not take higher orders into account 
because we are mainly interested in the qualitative behaviour of the freezing line for small 
densities p\. Denoting by F the Helmholtz free energy of the reference hard sphere system 
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(effective hard sphere diameter a), the total Helmholtz free energy F of the one- component 
system consisting of N\ big star polymers is approximated by 

^ = ^ + \fa J 90 (r)V pcTt (r)d'r. (19) 

In Eq. (119(1 above, g${r) denotes the radial distribution function of the reference hard-sphere 
system in the fluid phase and its angle-averaged counterpart in the solid phase, as defined 
in Ref. Q. 

We note here that more accurate methods for the treatment of potentials with a soft 
cores have also been proposed |6l| . but are not used here since we are only interested in 
the basic topology of the phase diagrams. For the free energy of the reference hard sphere 
system we used the equations of state of the Carnahan-Starling 



63( and Hall |64| for the 
solid and fluid phases, respectively. For the calculation of pair distribution functions we use 
the expressions of Henderson and Grundke j3| for the fluid and of Kincaid and Weis [f^ 
for the fcc-solid. We only considered the fcc-solid because this crystal structure appears 
at the fluid-solid transition in one component star-polymer solutions with an arm number 
h = 263 Q. 

B. Results 

In Fig. Owe show the resulting phase diagrams for size ratios q = 0.1 and q = 0.2 and 
functionalities of the smaller species fi = 16 and f'2 = 32, as obtained by the procedure 
described in the preceding subsection. The kinks for q = 0.1 at about p\<j\ = 0.8 and 
P2&2 = 0.1 and p\o\ = 0.6 and = 0.05, respectively, are an artifact of the choice 

(JTTj) and are associated with the sudden appearance of the Vp e rt(?")-term in Eq. (HHJ), once 
we cross over from the case d > d to the case d < d (cf. also Fig. EJ). Since we are 
primarily interested in the behavior for small densities p\ and the influence of the additives 
on crystallization, on the one hand, and on the possible existence of a spinodal line, on the 
other, a more sophisticated approach to the problem is at this stage not necessary. The 
liquid-solid coexistence region obtained by this approach is rather wide, due to the mapping 
on the effective hard-sphere system. Accurate calculations of the phase diagram of star 



polymers reveal that the coexistence region is much more narrow 



461 ] . yet the shape and 



evolution of the freezing lines as a function of p\ are not influenced by the width of the 
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FIG. 7: Phase diagrams for star polymer mixtures for different parameter combinations. The big 
star functionality is fixed at f\ = 263. The circles denote the phase boundaries as calculated by the 
hard-sphere mapping including the perturbation part V per t{r) and the lines are a guide to the eye. 
It can be seen that we only find a fluid-fluid demixing for q = 0.1. This binodal line is metastable 
with respect to the formation of the fcc-solid. The kink in the phase boundaries for q = 0.1 is an 
artifact and is caused by the method used to split the effective pair potential into reference and 
perturbation part. The symbol F stands for the fluid and the symbol S for the solid regions. 



density gap between the fluid and the solid phases. Finally, we note that we have shifted 
the freezing line to higher densities by an amount Ap±af = 0.06, in order to obtain the same 
density values for the crystallization as in the accurately known one component case [46]. 

For size ratio q = 0.1, the effective potential V e s(r) develops a strong attraction. This 
leads to a broadening of the coexistence area between the solid and fluid phase and eventually 
to a demixing binodal. This binodal is found, however, to be metastable with respect to the 
crystallization. The sharp kinks that show up in Figs. [7(a) and[7|[b) are artifacts of the way in 
which the effective hard sphere diameter a was determined and, in particular, of the fact that 
the attractive perturbation part l^,ert( r ) of the effective potential is absent in the treatment 
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for small reservoir densities below the kink and present above it. In reality, we expect the 
phase coexistence lines to 'turn around' smoothly, i.e., without the aforementioned artificial 
kink. However, the topology of the phase diagram and in particular the positive slope of the 
freezing lines and the subsequent broadening of the coexistence region into a 'gas-crystal' 
phase separation is not expected to be affected by these approximations. For larger values 
of q, a strong attraction does not emerge, so there is no fluid-fluid demixing. From Fig. 
13 it can be seen that less stars with /2 = 32 than with /a = 16 are needed to achieve 
similar effect. This is in agre ement with the properties of star-polymer-colloid mixtures 
which were investigated in [39]. The metastable binodal is closer to the stable region of the 
phase diagram for the smaller functionality of the depletant. The same trends were also 
observed in star-polymer-colloid mixtures 3^|, where stable binodals were only found for 
small depletant functionalities such as j'2 = 2 and f'2 = 6. For the star-polymer mixtures we 
consider here, the existence of a binodal is less likely than in star-polymer-colloid mixtures, 
because for star polymers the depletion force has to overcome the Yukawa-type repulsion 
between them before an effective attraction sets in. 

A striking effect is the melting of the crystal of the big stars upon addition of the small 
component, as can be seen from the positive slope of the freezing- and melting lines in Fig. 
|7| No effective attraction between the star-polymers is needed for this effect. As can be seen 
from Fig. El the effective hard sphere diameter a of the system decreases with increasing 
depletant density. Therefore, the effective packing fraction 77ns = (tt/6)pict 3 for the big 
stars gets smaller with increasing p\. Since this quantity determines the location of the 
freezing transition, the latter shifts to higher values of p\. The melting is therefore caused 
by the apparent shrinkage of the stars due to the depletion effects. The counterintuitive 
behavior of melting of a crystal although the overall polymer concentration is increased has 
its physical roots in the soft depletion mechanism and the associated reduction of the range 
of the repulsive potential between the big stars. 

In real experimental systems of star polymers, usually the formation of a glass is observed 
instead of crystallization. These glass transition lines are usually parallel to the freezing lines 
of the phase diagram, a property explicitly confirmed for the glass line of one-component 
star polymer solutions |50|. Therefore, we expect that the kinetic phase diagram of the 
mixture will have a topology running parallel to that of the equilibrium one that we traced 
in this study. We predict, therefore, that addition of small stars will bring about a melting 
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of the colloidal glass (or gel) formed by the large ones. This has been shown to be the case 
for mixtures of star polymers with linear homopolymer chains j^ j. 

V. COMPARISON BETWEEN ONE- AND TWO-COMPONENT DESCRIP- 
TIONS 

A. Chemical potentials 

By calculating the partial chemical potentials of the two-component system, we now map 
the reservoir representation (p 1; p 2 ) on the real physical system (pi, p 2 ), so as to make contact 
with experimental work, in which p\ has no direct relevance. The densities p\ and p 2 are 
linked by the condition that the chemical potential of the reservoir p\ has to be the same 
as the partial chemical potential in the system p 2 . 

The chemical potential in the reservoir of density p\ can easily be calculated if the fluid 
structure is known. The Helmholtz free energy density of the small stars in the reservoir, 
F/V = f(p 2 ), can be split into an ideal and an excess part: 

/(^) = /id(p r 2 ) + /cx(^), (20) 

where fid(z) = pln(^r 3 ) — z and r is an arbitrary length scale. The second derivative of the 
excess part is connected to the structure via the relation 66] 

f " M) = Pisikto-pi) ~ & (21) 

where S(k) is the structure factor of the small stars in the reservoir (one-component system). 
The free energy can then be calculated by two integrations, with the integration constants 
determined by the conditions f^ x (0) = and / ex (0) = 0. The chemical potential of the 
reservoir is given by 

/4 = f{pl)- (22) 

In order two calculate the partial chemical potential p 2 in the two-component system 
representation, we employ the so-called concentration structure factor 

Hx 2 2 S xl {k) + x\x 2 S 22 {k) - (x lX 2) 3/2 S 12 (k), (23) 
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FIG. 8: Mapping between the system- and reservoir densities p 2 and p\ of the small stars for 
f 2 = 32 and g = 0.1, and for various different values of the big star density p\. 

where Xi = pJ(pi + p-i)- Thermodynamic properties can then be calculated by using the 
equation 0, 0] 



lim S con (k) 



d 2 g(x 2 ,P,T) 



(24) 



where g(x 2 ,P,T) = G(x 2 ,N,P,T)/N is the Gibbs free energy per particle of the two- 
component mixture and P its total pressure. Eq. (|24|) can then be integrated as described 
in Ref. Bl to yield the sought-for Gibbs free energy per particle g(x 2 , P, T). Once the Gibbs 
free energy is known, the partial chemical potentials can be calculated using the equations 



g '(x 2 ) =H2-H\ 



(25) 



and 



g(x 2 ) -x 2 g'(x 2 ) = fa. (26) 

Reverting from the pair of variables (x 2 , P) back to (pi, p 2 ) for the mixture and using Eqs l22l 
l2Hl an d|2ni the mapping p T 2 ) — > (pi,p 2 ) can be carried out. Representative results are shown 
in Fig. |H| Clearly, the mapping depends parametrically on the big star density p\ in the 
mixture. For p x — one recovers p 2 = p\. In all cases, we obtain p 2 < p T s because all 
interactions are purely repulsive. The difference between the reservoir and system densities 
grows with increasing p x . 
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FIG. 9: Comparison of the fluid structure in the two- and effective one-component case. The 
densities are denoted in the diagrams. The structure of the mixture was calculated using the two- 
component RY closure, in the one component the usual RY closure was used. All plots are for 
q = 0.1 and ji = 32. In panels (a), (e) and (f) the solid and dashed lines fall on top of each other, 
so that the latter cannot be distinguished from the former. 

B. Structure 



We now consider the spatial correlations between the big stars in the fluid phase. We 
compare the correlation functions calculated in the two- component description with those 
arising from the one-component description using the effective interactions in the presence of 
the smaller species. The comparison has to be carried out for parameter combinations such 
that the reservoir- and system partial chemical potentials of the smaller species are equal 
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to one another. As the one-component description reduces the effective interaction to pair 
potentials, the comparison of the structural properties allows to estimate the magnitude of 
many body effects on the depletion. 

In Fig. IH1 we show the radial distribution functions and static structure factors of the 
big stars calculated with both different approaches. For the one-component case, the RY 
closure can be used because for the densities we consider here the effective interactions 
remain purely repulsive. The cases we consider here are for small densities p2 but a wide 
range of values of p\. The two component description includes many-body forces between 
the big stars which are caused by the smaller component. These effects are neglected in 
the effective one-component description j^| . In Fig. |§1 one can see excellent agreement of 
the results obtained by using these two different descriptions of the physical system. This 
additionally corroborates the validity of our approximation of the effective interaction. It 
can be also be concluded from the plots that the three-body forces are indeed much weaker 
than the pair interactions and can be safely neglected 



VI. SUMMARY AND CONCLUSIONS 



We have analyzed the structural and phase behavior of highly asymmetric mixtures of 
star polymers, with the asymmetry characterizing both their sizes and functionalities. The 
most striking phenomenon predicted by our investigations is the counter-intuitive melting 
of the colloidal crystal of the big stars upon addition of small ones. Though this finding ap- 
pears paradoxical at first sight, since addition of smaller stars increases the overall polymer 
concentration of the solution, its physical explanation can be traced to the effects of soft de- 
pletion. Whereas the depletants of hard spheres simply superimpose an effective attraction 
on a hard potential, when the big particles are themselves repulsive the depletion attraction 
is superimposed on a (soft) repulsion. In this way, a repulsive potential of reduced strength 
and/or range results and the effective, reduced repulsion is not any more sufficient to main- 
tain the stability of the crystal, which therefore melts. When the depletant concentration 
becomes sufficiently high, the attractive depletion force dominates over the soft repulsion, 
leading to a (possibly metastable) demixing transition between the two species. 

In real experimental systems of star polymers with high functionality, crystallization is 
hindered by the vitrification (gelation) transition and the large stars become structurally 
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arrested in a glassy state above the overlap concentration. The next step would be then to 
investigate the role and influence of smaller star additives on the glass transition of larger 
stars, a problem of significant importance for the control of the rheological behavior of soft 
matter through additives. Work along these lines is currently under way and the presentation 
of the results of this investigation will be the subject of a future publication. 
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